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Abstract. - We introduce and study two classes of Hubbard models with magnetic flux or with 
spin-orbit coupling, which have a flat lowest band separated from other bands by a nonzero gap. 
We study the Chern number of the flat bands, and find that it is zero for the first class but can be 
nontrivial in the second. We also prove that the introduction of on-site Coulomb repulsion leads 
to ferromagnetism in both the classes. 
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q Introduction. — Motivated by the recent discovery 
1 b f quantum spin Hall effect in band insulators with strong 
s'pin-orbit coupling (SOC) [IH3], the topological classifi- 
Gation of non-interacting electron systems has attracted a 
■^j- renewed interest [3H6] . The states of the systems are char- 
^sOacterized by the topological numbers linked to the pres- 
ence or absence of gapless edge modes carrying electronic 
^.or spin current. In integer quantum Hall systems [7], the 
I" first Chern number is directly connected to the quantized 
Hall conductance [8] , which is one of the most famous ex- 
amples of time reversal breaking insulators. On the other 
• • hand, in the recently found time reversal invariant insula- 
. £^tors, the states are classified by the Z2 topological num- 
ber. Since a topological number remains invariant as long 
$— 1 as the energy gap does not collapse, adiabatic transforma- 
^ tion from the original model to a flat-band model, where 
all the bands are dispersionless, provides a useful tool for 
the classification [2HS]- 



The flat-band models also play an important role in 
a completely different context, i.e., rigorous examples of 
fcrro- or fcrrimagnetism in the Hubbard model [9UT2]. In 
the models proposed by Mielke [10] and by Tasaki [11] , the 
on-site Coulomb interaction leads to ferromagnetic ground 
states when the lowest flat band is half-filled. A com- 
mon feature of these models is that the electron hoppings 
are frustrated and the lowest band is spanned by mod- 



erately localized eigenstates. Recently, similar localized 
states have been found in highly frustrated quantum mag- 
nets in strong magnetic fields jT3J[T3] and optical lattice 
models |15) . and offer a playground for studying nonper- 
turbative aspects of strongly correlated systems. 

These two subjects have developed separately, and the 
effect of electron correlation in topological insulators has 
not been studied intensively. Here we study the Hubbard 
model with flat bands in the presence of magnetic flux or 
SOC, which bridges the two subjects. In order to define 
the topological number, the gap between the flat band 
and other bands is required. The possibility of a gapped 
flat band is already a nontrivial issue since in many cases 
the band touching occurs at some points in momentum 
space [TS] and a uniform magnetic field destroys the flat- 
ness [17]. In this Letter, we propose two classes of tight 
binding models (TBM) which have a flat band separated 
from other bands by a nonzero gap [H] . The first class is a 
TBM on a line graph (e.g. checkerboard lattice) with non- 
uniform flux. For models in this class, we find that the flat 
band is always non-topological, i.e., the Chern number is 
always zero. We can study the effect of interaction rigor- 
ously and show that the ground states are ferromagnetic 
when the lowest band is half-filled for both magnetic-flux 
and SOC cases. The second class is a TBM embedded on 
a thin torus with a magnetic field perpendicular to the 
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Fig. 1: A graph consisting of 2 vertices and 5 edges, and cor- 
responding incidence matrix. The arrow indicates the sign of 
the phase factor e'° . 



plane. Surprisingly, all the bands become flat if a spe- 
cial condition is satisfied. We calculate the topological 
numbers and show that the topological flat band indeed 
exists. We also study the effect of electron correlation 
and rigorously show that the ground states are ferromag- 
netic when the lowest Landau level (LLL) is half-filled, 
i.e., v = 1/odd, without using the LLL projection |19j . 

Non-topological flat band. — We start from the 
first class. Let G = (V, E) be a graph (lattice), where V is 
the set of vertices (sites) and E the set of edges (bonds) . 
We assume that G is twofold connected [20] . We define the 
incidence matrix B = (B ve ) ve v.eeE by assigning a nonzero 
complex number to B ve when the vertex v is incident to 
the edge e, and by setting B ve = otherwise (see Fig. Q] 
for example). Define the line graph L(G) = (Vl,-El) 
of G as usual, by first regarding (the midpoint of) each 
edge in E as a vertex in Vl (thus Vl = E), and then 
connecting any pair of vertices in Vl (by an edge in E^) 
when the corresponding pair of edges in E share a com- 
mon vertex. Fig. [5] (a) shows the square lattice and its 
line graph, the checkerboard lattice (ignore the differences 
in the bonds for the moment). Let us consider TBMs on 
G and L(G) with hopping matrices (single-particle Hamil- 
tonians) T = BB^ and Tl = B^B, respectively. It is eas- 
ily shown that (Tl) all the eigenvalues of T and Tl are 
nonncgative, (T2) nonzero eigenvalues of T and Tl are 
identical, and (T3) Tl has at least (\E\ — |V|) zero-energy 
eigenstates [5T]. (Tl) follows if one notes that both T and 
Tl are positive semidefinite. (T3) is an immediate conse- 
quence of the fact that B is a \V\ x \E\ matrix. To see (T2), 
let ip be an eigenvector of T with a nonzero eigenvalue a, 
i.e, BBt<p = cap. Multiplying from the left, one finds 
Tl^> = cap with a nonzero vector ip = B^ip. To complete 
the proof we only need to repeat the same argument with 
T and Tl switched [25]. For periodic systems, these zero 
energy eigenstates of Tl form the lowest flat band. 

Let us now show that the flat band in the TBM on L(G) 
may be gapped but is not topological. To be concrete we 
let G be the square lattice [53], and consider models with 
a constant flux (j> per plaquette. For each edge (vv') E E, 
define <j> vv i = —(f) v ' V £ M so that <j> VlV2 + <j> V2V3 + 4>v 3 v 4 + 
4>vim = 4> (mod 1) for any plaquette (v\V2VzVa) oriented 
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Fig. 2: (a) Left: Tight-binding model on a square lattice 
with magnetic flux cj) through each plaquette. Solid and broken 
edges belong to E' and E" , respectively. The rotating arrow 
indicates the orientation of vertices. Right: Line graph of the 
left graph. The flux <f> goes through diamond-shaped plaque- 
ttes. The hopping amplitudes |(TL(a?)) ee '| of solid, dotted, and 
broken edges are f, y/x, and x (0 < x < 1), respectively, (b) 
Honeycomb lattice and its line graph (Kagome lattice). Nota- 
tions are the same as in (a). 



in the counterclockwise direction. By setting B ve = 
exp[m(j) vv i ] if e = (vv'} , we get a TBM on G with a uniform 
flux, and the corresponding TBM on the line graph L(G} 
with a uniform flux through the diamond plaquettes (see 
Fig- [11(b), still ignoring the differences in the bonds) |24j . 

It is useful for our proof to introduce interpolating 
TBMs with an extra parameter < x < 1. Consider a 
disjoint decomposition E = E' U E" as in Fig. [T|b), and 
redefine B ve as y/x exp[wi(p vv i] only if e = (vv'} € E" . 
We denote by T(x) and T\,(x) the corresponding hop- 
ping matrices. Note that T(l) and Tl(1) are the same 
as the original T and Tl, respectively, and both T(0) and 
Tl(0) describe TBMs which decouple into local pieces. 
From the definition, one finds for any < x < 1 that 
T(x) > T(0) > e(4>) where e(4>) is the lowest eigenvalue of 
T(0) [5S]. Suppose that <j> is nonintegral. Since an explicit 
calculation shows e(</>) > [55], the above (T2) implies 
that Tl(x) has a nonzero gap above the zero eigenvalue. 
Recalling that the lowest flat band of Tl is gapless for 
<j> = 0, we see that the above gap in Tl originates from 
the flux. 

Now we investigate the Chern number defined by impos- 
ing twisted boundary conditions as in [571 [55] of the flat 
band in the TBM on the line graph L(G). Suppose that 
<f> is nonintegral. Since the Chern number is invariant as 
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Fig. 3: A portion of the Kagome ladder. Unfilled circles are 
vertices of the square ladder. The shaded region depicts the 
localized state /' and its amplitudes are indicated with r — 
e l6 / 2 . Fluxes through plaquette and triangle are —0/(2n) and 
0, respectively. 



long as the energy gap does not close, one can evaluate it 
in the TBM with Tl(x) for any x. But since a decoupled 
system is insensitive to a twist in the boundary conditions, 
the Chern number of the flat band in Tl(0) is clearly zero. 
This proves that the Chern number of the flat band in the 
original model with Tl is zero. 

Although the fiat band is non-topological in this class 
of models, we can rigorously study the effect of electron 
correlation. We present two specific examples in the 
following. 

i) Kagome ladder. — The first example is the Hub- 
bard model on the Kagome ladder, the line graph of the 
square ladder, shown in Fig. [3J The sites in the mid- 
chain are labeled by e = nd with n = 0, 1, . . . , 2L — 1 as 
shown in Fig. [3J We denote vertices of the square ladder 
by (e, /) with I = ±1. We let the incidence matrix ele- 
ment -B( e ,z) e / be 1 if e! = e or e — 11/2, e l8 ^ 2 if e' = e + lu\, 
and otherwise. Here 8 is the parameter determining the 
flux, and Ui and are the vectors indicated in Fig. |3J 
The tight-binding Hamiltonian of our model is then given 

by 

<t=T4 o-=t,4- 

l=±l 

(1) 

where e in the right-hand side is summed over the 2L sites 
in the mid-chain, and the a-opcrators arc defined as 

a (e,/),cr ^ ^ -B(eJ) e'^e' ,a ~ ^e,cr H~ e ^ Ce-\-lvi ,<r H~ c e — jCr . 
e' 

(2) 

We impose periodic boundary conditions in the chain di- 
rection (see Fig. |3j. By a straightforward calculation, one 
finds that 

/,U = -(e 4e + l) C L 

+ V e L e+lv l ,a + C e+lv 2 ,<J e C e+lS,a) \°) 

l=±l 

anticommutes with the a-operators. Therefore, the single- 
electron zero-energy states are given by /^J^o), where 



|<J>o) is the vacuum state. The collection of these states 
forms a complete basis for the fiat band. 

We shall consider the case where the flat band is half- 
filled. In the non-interacting case, the 2L-electron ground 
states are highly degenerate and exhibit paramagnetism. 
Let us add the Hamiltonian the standard on-site Coulomb 
repulsion 

H v = u ^2 n e . t n e ^, (4) 
eev L 

where n e (T — c\ 5 c C II is the number operator at a site 
e of the Kagome ladder. Then the degeneracy is lifted 
and only the ferromagnetic states remain as the ground 
states, which are the zero energy eigenstate of both Hkl 
and Hjj. To prove this claim, we introduce new fermion 
operators 4,<r = e^irU ~ /l+(_i)n >CT , and follow the 
standard strategy [25]. The states dj^l^o) form an- 
other basis for the flat band. By representing a ground 
state |\E') in terms of the d-operators, we can firstly show 
that the zero-energy conditions for the on-site repulsion, 
Ce'^Ce'^l^) = 0, with sites e' = nS + {—l) n u\ forbid the 
double occupancy of d-states. Then the same conditions 
with sites in the mid-chain imply that the ground states 

must be (n^o 1 ^n,t) l^o) an d its SU(2) rotations. We 
note that even when the flat band is less than half-filled, 
ferromagnetic states are ground states but are not unique. 

It is also possible to consider the non-topological flat 
band model associated with SOC. For the Kagome ladder, 
the corresponding tight-binding Hamiltonian i?KL * s P r0 ~ 
duced by replacing e l6 / 2 in the definition of a-operators 
with e" e / 2 in Eq. ([T]). This model has spin-dependent 
complex hoppings. Since for = 7r, the model is reduced 
to the case of spin-independent hopping by a local gauge 
transformation, we restrict ourselves to 6 7^ it. The single- 
electron zero-energy states of H§^ c are given by ct^ a \$o), 
where a is defined as (y n a with 

/L = -(^ 9 + i)4,. 

+ V (V^/V 4- c x - P iCT< VM ) (5) 

+ V C e+lu u a + C e+lu 2 ,a e C e+lS,a) \° ) 

!=±1 

in place of /t a . In the 2L-electron case, the ground 
states of the Hubbard Hamiltonian H = -ff^L + are 
given by frinio 1 ^n,aj l^o) with a =t,4- In contrast to 
the spin-independent case, the SU(2) spin degeneracy is 
lifted and there are only two ground states. 

ii) Two-dimensional checkerboard lattice. The 

second is the Hubbard model on the checkerboard lattice, 
the line graph of the square lattice whose vertex set is 
given by V = [0,L— I] 2 n Z 2 with an odd positive integer 
L. The periodic boundary conditions are imposed in both 
directions. We label an element of the edge set E of the 
square lattice by its mid-point position. Let fix = (1/2, 0) 
and H2 = (0, 1/2). For later use we decompose E as E = 
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Fig. 4: Checkerboard lattice (Line graph of the square lat- 
tice). Filled (unfilled) circles represent sites corresponding to 
the edges in E2 (£a). The shaded region depicts the localized 
state d ei(T in Eq. 



E1UE2, where E\ = {e = v+fii\v E V}. Let the incidence 
matrix element B ve be 1 if e = v ± H2 or e = v — fix, 
e i2(v/j, 2 )0 \f e = v + fii, and otherwise. Here 9 is again 
the parameter determining the flux (f> = —9/(2tt). The 
tight-binding Hamiltonian of this model is then given by 
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cl 



E 



' C e,ir C e 



CT =t4 

e.e'GVL 



E 

0=t4 



(0) 



where the a-operators are defined as 



e i2(vfj, 2 )e c 



— fj, 1 , a 



(7) 



For each e S E\ we define 



i(2e»jn 2 ) 



L-l 
i=0 



(-i)'4- M - 



2//X2 :C 



L-l 

E( 

Z=0 



-lVc 



e+ fj.+2l jj, 2 



OS) 



where /x = /^i + /^2- These states are extended in one 
direction and localized in perpendicular one as shown in 
Fig. I?] It is easy to see that the d-opcrators are anti- 
commute with the a-operators, and therefore the single- 
electron zero-energy states of Hql are given by d\ a \<&o). 
The collection of these states forms a complete basis for 
the flat band. 

As in the case of the Kagome ladder, the fully polar- 
ized states, ^IleGBi °C t) l^o)' and its SU(2) rotations 
are the unique ground states of the Hubbard Hamilto- 
nian H = Hql + Hjj when the electron number is L 2 , i.e., 



the flat band is half-filled. Note that these ground states 
are simultaneous eigenstates of both ifcL and Hjj with 
zero-energy. This claim can be proved by following the 
same strategy. Representing a ground state in terms 
of the d-opcrators and noting that e in E\ supports only 
d ei(T , we can firstly show that the zero-energy conditions 
for the on-site repulsion, c e _^c e ^\ty) = 0, with sites e£ J?i 
forbid the double occupancy of d-states. Then, examing 
the same conditions with sites e G E2 we arrive at the 
conclusion. 

Topological flat band. — We turn to the second 
class. We consider the square lattice TBM embedded on 
a torus with a magnetic field perpendicular to the plane. 
Such a problem is known as the Hofstadtcr problem |30j 
and has been extensively studied [31ti33] . As we will show, 
all the bands become flat if the flux per plaquette and the 
number of sites along (1,1) direction satisfy certain con- 
ditions. We shall use a notation as close to those in [35] 
as possible. The tight-binding Hamiltonian is given by 

-^hop T x 



with 



Tr 



EE 

CT=f4 m ' n 

EE' 

er="|\l m , n 



(10) 



where (m, n) denote the vertices of the square lattice, and 
®m n = ( m + n ) 7T 4 > and 9m n = — ( m + n + l) 71 "^- The flux 
per plaquette is <j) = P/Q with mutually prime P and Q. 
Periodic boundary conditions are imposed both in (1,1) 
and (1,-1) directions. From the Bloch theorem, we can 
assume the single-electron state to be of the form 



^* m , n (p +) p)cJ mra)jiT |$o>, (11) 



l (p+,p) - e^( m+ ") + ^™-"V™ + „b + ,p),(12) 



where ip tc+2 Q(p+,p) = ^ k (p +1 p), (k = 0, 1, 2Q- 1). We 
now consider the thin torus case where the system size is 
finite along one of the cycles of the torus (sec Fig. EJa)). 
Suppose that the periodic boundary conditions in (1,1) 
and (1,-1) directions are given by 



^m+Q,n+Q(p+,p) = $ m ,n(P+>P) 



S S> m +L,n-L(p+,P) 



®m,n(P+,P)- 



(13) 
(14) 



Then they yield p+ = nl/Q (mod 7r/Q) with I S Z, and 
p = kir/L with k = 0, 1, L — 1. Therefore, if Q is even, 
p + = tt/2 (mod tt/Q) and it automatically satisfies the 
mid-band condition. This condition is closely related to the 
symmetry described by the quantum group U q (sl2) [311132] . 

Henceforth we shall focus on the case of even Q and 
show that all the bands are flat, namely, the single-electron 
energy e{p) is independent of p. The Schrodingcr equation 
for |$ cr (7r/2,p)) is written as 



i(q 



i+i _ q -(l+l) 



)u l+ i+i(q-q )ui-i=c{p)u h (15) 
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Fig. 5: a) Thin torus. Independent cycles in (1, 1) and (1, —1) 
directions are indicated by the broken lines, b) Lattice struc- 
ture for Q = 2 and L = 3. (m, n) denotes the label for the 
vertex. Periodic boundary is indicated by the dashed line, c) 
Table for the sum of the Chern numbers of the lowest flat 
bands. n c and v are the electron density and the filling factor, 
respectively. The unfilled circles indicate that the Chern num- 
ber is undefined since there is a gap closing caused by a twist 
in the boundary conditions, (d) Single-electron band structure 
for P = 1 and Q — 4. Each band is doubly degenerate (per 
spin). The sum of the Chern numbers of the lowest bands is 2. 



where q = e mP /® , and ui are defined through the unitary 
transformation 



1 2Q-1 
ip k {n/2,p) = -j= <1 



Ik e- Up ui. 



(16) 



i=o 



It is now obvious that e(p) does not depend on p because 
there is no p-dependence in the LHS of Eq. (TT5"]) . The 
single-electron energy for P = 1, Q = 4 is shown in Fig. 
O (d) . Each band is doubly degenerate for the reason dis- 
cussed below. Interestingly the lowest flat band in this 
construction can be topological. In Fig. \E\ (c), we list the 
Chern numbers for several 4> = P/Q computed numeri- 
cally using the method of [271125] . 

We next show that the spatially localized state along 
(1,-1) direction can be constructed from the solution of 
Eq. (|15p . Due to the fact q® — q~® = 0, we can take 

{^}f=o _1 t0 bc of the form: ui = ui 6 R if < K 
Q — 1 and zero otherwise. This solution is degenerate 
with the other one: ui = (—lyvi-Q if Q < I < 2Q — 1 and 
zero otherwise [33]. The vector {vz}^ 1 is normalized as 

Y^Lo 1 v t = 1- From those solutions, a localized Wannier 
state extending from m — n = j to j + Q — 1 can be 
constructed as 



4.1^) = 



Q-12Q- 



l*o>, (17) 



where c, . = c] k+i k _.. 

the same parity and otherwise. One can easily show 
that {dj^,d\ >tr ,} = 5ji6 aa - (0 < j,l < 2L - 1). From the 
Perron-Frobenius theorem, one also finds that the lowest 
eigenvalue for Eq. ([15]) is two-fold degenerate. This im- 
plies that the lowest energy of Hi wp is 2i-fold degenerate. 
Therefore, d-states form a complete basis for the lowest 
flat bands. 

Let us now study the effect of electron correlation within 
the Hubbard model. We define the Hubbard Hamiltonian 
H by 

H = H hop + E n (™»»).t n (™,n),-l (1®) 

with U > 0. If the total number of electrons N e is 2L, 
the ferromagnetic state constructed from the lowest flat 
band of -Hhop, |^) = E[j=o 1 dj -f-l^o), is a ground state of 
H. To go further and show that all the ground states are 
ferromagnetic, we make use of theorem due to Mielke [35] . 
The theorem asserts that if the single-particle density ma- 
trix constructed from as 



1 



Pr,r> = ly-(*|c{ roi „ ))t C( m% „'), t |*) 



(19) 



is irreducible, I'S) is the unique ground state of H with 
N c = 2L up to the trivial degeneracy from the SU(2) 
symmetry [36] . Here, r and r' denote (m,n) and (m',n'), 
respectively. The diagonal matrix element, the electron 
density, is uniform and obtained as p r _ r = 1/(2QL). To 
show that p r r i is irreducible, it is sufficient to show that 
any matrix element corresponding to a pair of nearest 
neighbor sites is nonzero. Moreover, p r r i is Hcrmitian, 
we have only to study the case of m' — n' = m — n + 1, 
vnl + n' = m + n ± 1. In this case, an explicit calculation 
gives 



^ i±V ( ' +1 Wi. 



2QL 



1=0 



(20) 



1=0 k=0 



Recalling that {wzj^Q 1 satisfy Eq. (TT5"]) . we find 
f^ c EE=o 2 * ±1 9 ± ^ +1 ^iWi+i] = ei/4, where ei is the single- 
particle energy of the lowest flat band. For even Q, e± ^ 
can be shown by deriving the characteristic equation for 
e(p) from Eq. (TT5]) and hence p r y ^ 0. Therefore, via 
Mielke's theorem, the ferromagnetic state |^) is the unique 
ground state of H up to the spin degeneracy. This ground 
state corresponds to the quantum Hall ferromagnet for the 
filling factor v = 1/P with odd P. So far, we have stud- 
ied the square lattice model on the thin torus. However, 
our argument can be generalized to other lattice geome- 
try such as the honeycomb lattice which is relevant to the 
quantum Hall ferromagnctism in graphene [37] . 

Conclusion. — To conclude, we have studied two 
classes of Hubbard models with a lowest flat band sep- 
arated from the other bands by nonvanishing gap origi- 
nating from the magnetic flux or spin-orbit coupling. In 
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the first class, we have shown a systematic way to con- 
struct tight binding models with gapped flat bands using 
line graphs. We have proved that the Chern number of 
the flat bands is zero in this class of tight binding models. 
We have also studied the effect of the on-site Coulomb in- 
teraction for two particular cases: i) Kagome ladder and 
ii) Two-dimensional checkerboard lattice, and have rig- 
orously shown that the ground states are ferromagnetic 
when the lowest flat band is half-filled. 

In the second class, we found the construction of the 
tight binding models embedded on a thin torus, in which 
all the bands are flat. Each flat-band manifold is spanned 
by the states localized in one direction while dclocalized in 
the other. This is reminiscent to the LLL wave functions 
on a torus. We have numerically calculated the Chern 
number of the lowest band and found that it can be non- 
trivial. The lowest flat bands also allow us to study the ef- 
fect of the on-site Coulomb interaction nonperturbatively. 
Applying the theorem of Mielke, we have shown that the 
ground states are ferromagnetic when the lowest flat band 
is half-filled. Although our model only reproduces inte- 
ger quantum Hall systems, it would be interesting to ex- 
plore lattice realizations of fractional quantum Hall sys- 
tems where the next nearest neighbor interaction is prob- 
ably important |38j . 
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